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We rigorously apply the sum rules to the sum-over-states expression to calculate the fundamental 
limits of the dispersion of the two-photon absorption cross-section. A comparison of the theory 
with the data suggests that the truncated sum rules in the three-level model give a reasonable 
fundamental limit. Furthermore, we posit that the two photon absorption cross-section near the 
limit must have only three dominant states, so by default, the three-level model is appropriate. This 
ansatz is supported by a rigorous analytical calculation that the resonant term gets smaller as more 
states are added. We also find that the contributions of the non-explicitly resonant terms can not be 
neglected when analyzing real molecules with many excited states, even near resonance. However, 
puzzling as it may be, extrapolating an off-resonant result to resonance using only the resonant term 
of the three- level model is shown to be consistent with the exact result. In addition, the off-resonant 
approximation is shown to scale logarithmically when compared with the full three-level model. This 
scaling can be used to simplify the analysis of measurements. We find that existing molecules are still 
far from the fundamental limit; so, there is room for improvement. But, reaching the fundamental 
limit would require precise control of the energy-level spacing, independently of the transition dipole 
moments - a task that does not appear possible using today's synthetic approaches. So, we present 
alternative methods that can still lead to substantial improvements which only require the control of 
the transition moment to the first excited state. While it is best to normalize measured two photon 
absorption cross-sections to the fundamental limits when comparing molecules, we show that simply 
dividing by the square of the number of electrons per molecule yields a good metric for comparison. 



I. INTRODUCTION 

Two-photon absorption is the fundamental pro- 
cess behind most resonant applications such as 
three-dimensional photolithography, Q] photodynamic 
therapy, [3 high-density storage [3 and optical power 
limiters.Q All such applications require a large resonant 
two photon absorption (TPA) cross-section. The devel- 
opment of a deeper understanding of the TPA cross sec- 
tion could be used as a guide to optimizing the nonlinear 
response of molecules. Our approach is to use sum rules, 
which relate transition dipole moments and energies to 
each other, to learn how the various terms in the sum- 
over-states expression interfere with each other and how 
molecules can be designed for constructive interference. 

The maximum fundamental limits of the largest di- 
agonal tensor component of the off-resonance first and 
second hyperpolarizabilities have been calculated using 
the generalized Thomas-Kuhn sum rules. 0, 0> Ex- 
trapolating this off-resonance result of the second hy- 
perpolarizability made it possible to estimate the maxi- 
mum values of the resonant two-photon absorption cross- 
section. [j| Using this process, it is found that most or- 
ganic molecules show a TPA cross section that falls 
well below the estimated limit (by a factor of at least 
1CP 3//2 ), suggesting that the paradigm for making organic 



molecules for TPA applications may not be optimized. 
This factor appears to be universally observed for single 
molecules in all nonlinear processes studied (such as sec- 
ond harmonic generation, hyper rayleigh scattering. [TEj 
and two photon absorption) ;[9j but can be breached if 
the molecules are made to interact. I n anv case, how- 
ever, the fundamental limit has not been breached, nor 
is it expected that it would be possible to do so since the 
sum rules are as basic as the Schrodinger equation from 
which it derives. 

In the present work, we apply the sum rules directly to 
the full dispersion of the TPA cross-section without ex- 
trapolation and without neglecting any terms. As such, 
these results are exact, in analogy with past calculations 
of the fundamental limits of the off-resonant hyperpo- 
larizabilities. We use our approach to analyze a series 
of molecular measurements which shows that the fun- 
damental limits are somewhat lower than predicted us- 
ing extrapolation from the off-resonant result; but, the 
molecules analyzed are still well below the fundamen- 
tal limits. More importantly, we show that Two-Photon 
Fluorescence measurement, used in getting TPA cross- 
sections, may be inaccurate if all states are not included 
in the analysis. 
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II. THEORY 

The TPA cross-section S (in units of cm 4 s) is related 
to the imaginary part of the diagonal component of the 
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second hyperpolarizability, 7/, through: [b"| 
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where 7/ = 7771(7) = ^(7a;xxa(-w; w, w, -a;)), 77 is the 
refractive index of the bulk material, u) the frequency of 
the incident light, c the speed of light, and the brackets 
indicate the average over all-possible orientations of the 
molecule. For a one-dimensional molecule (7) =7/5 and 
for a spherical molecule (7) = 7. The effects of the local 
fields are included in the dressed hyperpolarizability, 7J . 

Two photon absorption is inherently a resonant pro- 
cess because it quantifies the strength of absorption of 
two photons - each with energy ftw - from the ground 
state of energy E to the second excited state with en- 
ergy E2, where E213 = E% — Eq = 2Hui. Consequently, 
the exact expression for the diagonal component of the 
second hy perp olarizability given by the sum-over-states 
expression .ll) is often approximated by the explicitly res- 
onant term in the sum, "f res : 
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where E n o is the transition energy between the excited 
state 77 and the ground state, \i nm is the transition dipole 
moment between states 77 and 777, and r n0 is the damping 
factor (inverse radiative lifetime) between the state 77 and 
the ground state. When the approximation 7 w j res is 
made, the contributions of all the non-explicitly resonant 
terms to 7 are ignored. For such an approximation to be 
valid, 7 res should dominate the response in the domain 
of interest for TPA processes. 

Two different regimes can be distinguished in TPA pro- 
cesses: the single resonance regime and the double res- 
onance regime. In both cases E20 rs 2hco] but, in the 
single resonance regime, the photon energy fko does not 
match the transition energy Eiq (i.e. \Eiq — hu\ ^> Tio ) 
while in the double resonance regime Eiq ~ hw. For 7res 
to dominate the second hyperpolarizability 7, it is neces- 
sary to work in the single-resonance regime to ensure that 
j res is the only "resonant" term in the sum-over-states 
expression. 

Studies of the first and second hyperpolarizabilitics 
in the off-resonance regime clearly show that individual 
terms in the sum-over-states (SOS) expression cannot 
be independently adjusted since the terms are related 
to each other through the sum rules. [EEl In fact, in the 
SOS expression for 7 there are additional two-photon res- 
onant terms similar to Eq. [3 that contribute to 7. Even 
if certain terms appear to dominate the response, all the 
contributions must be included since the contribution of 
just one term might not be representative of the full re- 
sult - especially if there are cancellations between large 
terms. Therefore, our calculations of the maximum limit 
of the TPA cross-section includes all contributions to 7. 
We also evaluate the relevance of the different terms in 
the SOS expression with the aim of determining whether 



or not Eq. [3 is in general a good approximation to 7/; 
an approximation that is used in determining the TPA 
cross-section (and /112) with TPA fluorescence measure- 
ments. 

The exact sum-over-states expression for 
lxxxx{— w; u), u), —lu) Hj can be split into two types 
of contributions as follows: 

Ixxxx (-uj;uj,uj, -w) = Si+S 2 +S 3 +S4,+Ti+T 2 +T 3 +T i , 
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Here, /123 denotes the average over all distinct permu- 
tations of w\, u>2 and u> 3 , where ojj = lu, u)% = lu and 
lu 3 = —lu. The prime in the sum indicates that the sum 
is over the excited states, Q n o is a complex quantity de- 
fined as: 
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Note that the explicitly two-photon single-resonant terms 
(when E20 — 2Hlu) are all contained in S\ and S2 (Eqs. 
0] and [5J. Ti and T 2 (Eqs. [H] and EJ) contain terms that 
will be resonant in the double-resonance regime (when 
_Bio = huj and E20 = 2fku). 

It is important to keep sight of the approximations in- 
herent in the SOS expression, which is at the heart of 
all our calculations. It assumes that the light emitted 
by a molecule is electric dipole in nature, so magnetic 
transitions as well as higher-order electric terms are ig- 
nored. In addition, since the SOS expression is derived 
from perturbation theory, the energy of interaction be- 
tween the light and the molecule must be small compared 
with the binding energy. As such, the light can not excite 
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the molecule into a real high-lying state. Finally, we fo- 
cus on purely electronic transitions, so vibronics are not 
considered. However, it has been shown that neglect of 
nuclear motion has only a small affect on the limits of 
the calculated nonlinear-optical response.. 15f 

We consider a general three-level model and truncate 
the Thomas-Kuhn sum rules to the first three energy 
levels to obtain relationships between the different pa- 
rameters that determine 7. The procedure is similar to 
the one followed to obtain the maximum value of the off- 
resonance first and second hyperpolarizabilities.0,0] Un- 
like this past work, the present work includes the effects 
of dispersion, so the widths of each state are included 
by using an imaginary part of the frequency, as shown 
in Equations 0] through ^] Hence, the expressions are 
mathematically more algebraically messy. We therefore 
use MATHEMATICA® to evaluate the results. 

The full set of self-consistent three-level truncated 
Thomas-Kuhn sum rules yield relationships between 
transition dipole moments and energy differences as fol- 
lows: 
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where A/j, m „ = (i mm -/i nn , e is the charge of an electron, 
N is the number of electrons in the molecule and m is 
the mass of an electron. 

A direct consequence of the first sum rule given by Eq. 
1141 is that the values of \/J,io\ are constrained: 
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Therefore, as in the off-resonance analysis, we define the 
dimcnsionlcss quantity, X: 



X 

where -1 < X < 1. 
quantity E as: 



IM10I 



I,, MAX 

IM10 



< 1, 



(19) 



We also define the dimcnsionlcss 
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with < E < 1. We note that this explicitly assumes 
that E10 < E20 and that -E20 is the two-photon state. 
These criteria are obeyed for all of the molecules ana- 
lyzed in this paper. Note that systems such as the longer 
polyenes have lower-energy two-photon states. As such, 
a separate calculation in the spirit of the one that fol- 
lows would need to be developed. While we do not do 
such a calculation in this paper since it is not relevant to 
our studies, interesting results are possible. Secondly, we 



note that if the first two excited states in an asymmetric 
molecule become degenerate, then the next excited state 
will need to be considered and the energy ratio redefined, 
i.e. E — > E10/ E 30 - 

Our objective is to express 7/ in terms of X and E 
and determine the fundamental limits of the TPA cross- 
section, S(ui), by finding the values of X and E that 
maximize 5. It is important to emphasize that the calcu- 
lation of the fundamental limits makes only one assump- 
tion: the forces between the charges are described by a 
potential energy that depends exclusively on the coordi- 
nates of the electrons and the nuclei. The tools that we 
use to calculate the fundamental limits, the generalized 
Thomas-Kuhn sum rules, are fundamental identities de- 
rived from the Schodinger equation, and are applicable 
to all electronic excitations, regardless of the particular 
details of the molecule. However, we approximate the 
molecule as a three-level system to model the response 
and truncate the sum rules to three states, an approach 
that has been successful at determining the fundamental 
limits in the off-resonance regime. In fact, recent re- 
search suggests that molecules with more excited states 
have a smaller first hyperpolarizability than the three- 
level system. ^5| As such, the limits we calculate are truly 
an upper limit. 



A. Contribution of the explicitly resonant 
two-photon terms 

We first evaluate the contributions of the terms Im(Si) 
and Im{S2) to 77, since they contain all the two-photon 
resonant contributions (when E20 = 2hu>). Using the 
three-level-truncated Thomas-Kuhn sum rules, we are 
able to parameterize Im{S\) and Im(S2) in terms of £10, 
-E201 Xt Tin, and I^o- Since these expressions are con- 
strained by the three- level Thomas-Kuhn sum rules, they 
are labelled Sf K3 . We assume at all times that the two- 
photon resonance condition is obeyed, that is -E20 = 2Huj; 
and we vary the other parameters. 

We find that both Im{S^ K3 ) and Im{S% K3 ) vanish 
when X = and reach their maximum value when X = 
1. Also, as Tio and I^o increase as a fraction of Eiq, 
the values of Im(Sf K3 ) and Im(S'J K3 ) approach zero. 
For our calculations we consider only one electron (N=l) 
and we use typical values of the parameters for organic 
molecules: 



X = 0.5, 
Tio = r 20 = O.leV, 
and E 10 = 2eV. 



(21) 
(22) 
(23) 



Fig. UJshows a plot of Im{Sl K3 ) and Im{S' 2 TK3 ) indi- 
vidually as a function of the energy ratio E — E\q/E2q- 
Also plotted is the behavior of their sum, Im(Sf K3 ) + 
Im(S^ K3 ), and Im{ lres ) (Eq. 0. Both Im{S? K3 ) and 
Im(S2 K3 ) are peaked at E = 1/2, and fall quickly to 
zero as E — > 1 and E — > 0. As we later discuss, this 
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FIG. 1: Im(ST K3 ), Im(S% K3 ), Im(S 1 + S 2 ) and Jm( 7res ) as 
a function of the energy ratio, E, with X = 0.5, Fio = T2o = 
O.leV, E 10 = 2eV, and N = 1. 



FIG. 2: Im(ST K3 ), Im(S% K3 ), Im{T? K3 ) and Jm(T 2 TK3 ) as 
a function of the energy ratio, E, with X = 0.5, Fio = T2o = 
O.leV, E 10 = 2eV, and N = 1. 



is the general trend for all the terms that contribute to 
7/, resulting in an imaginary third-order susceptibility 
that peaks around E = 1/2. As such, to maximize the 
TPA cross-section, the two lowest excited state energies 
of a molecule must be equally spaced to allow for the 
double-resonance condition (E 2 o = 2£io = 2^o;).[l^ 

While individually Im{Sj K3 ) and /m(S'^ if3 ) are 
peaked functions at E = 1/2, they have opposite signs, 
so these two terms partially cancel. More importantly, 
although Im(Sf + S^ 3 ) contains the term 7m( 7r . es ), 
it is obvious from the spectrum predicted by Im(Sf K3 + 
S% ) that it is quantitatively and qualitatively differ- 
ent from Im^res). Furthermore, Jm(7 res ) diverges as 
E — ¥ 1, while Im(Sf K3 ) and Jm(S'J K3 ) approach zero. 
When all terms are included, the divergences from indi- 
vidual terms in Im(Sf K3 ) and Im(S2 K3 ) cancel in the 
same way as it does in the off- resonance calculations. 0@ 



B. Full three-level calculation of 7/ 
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FIG. 3: Im(SZ K3 ), Im(Sj K3 ), /m(T 3 TK3 ) and 7m(T 4 TK3 ) as 
a function of the energy ratio, E, with X = 0.5, Fio = T20 = 
O.leV, E10 = 2eV, and N = 1. 



The analysis of the sum Im(Si K3 ) + Im(S2 K3 ) shows 
that to get a substantial TPA cross-section the system 
must be close to the double-resonant condition (when 
E20 = 2£ , 1 o = 2huj). In this regime, contributions of 
terms such as Im(T^ K3 ) and 7to(T|" K3 ) can play an es- 
sential role, since they explicitly contain one-photon res- 
onant terms (E\a = fko). Therefore, we now study the 
behavior, individually, of all the terms that contribute 
to the second hyperpolarizability as constrained by the 
three-level truncated sum rules, r yJ K3 - We plot each term 
(defined by Eqs. 1411 Q as a function of E in Figs. [5] and 
13 For simplicity we have used the values defined in Eqs. 
I21ll22l andl23lsince these are the typical values for organic 
molecules that absorb light in the visible. 

Before proceeding, we must take a small detour to deal 



with the fact that on double resonance, we must con- 
sider the other excited states that may contribute to the 
SOS expression for 7/. This is certainly a crucial issue 
when modelling real molecules, which have many excited 
states beyond the first three. However, we must also re- 
call that our analysis focuses on understanding the fun- 
damental limits. Thus, our calculations are based on the 
ansatz that the hyperpolarizability, of any order, will be 
maximized for any system when all of the transitions are 
concentrated into three states. As such, when we are 
dealing with ideal molecules that are near the fundamen- 
tal limit, the problems of higher lying states contributing 
to 7/ is moot. The next section addresses this issue for 
the explicitly-resonant term. 
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FIG. 4: yf Ki as a function of E and X, with I\ = T 2 o = 
O.leV, E 10 = 2eV, and N = 1. 



While the contributions of Im(Sj K3 ), Im(Sj K3 ), 
Im(T^ K3 ) and 7m(T 4 Tif3 ) are negligible m comparison 
with Im(Si K3 ) or Im(S2 K3 ), it is obvious from the 
plots that the contributions of Im(T^ K3 ) and Im(T^ K3 ) 
are comparable to the contributions of Im(S-[ K3 ) and 
Im{S^ K3 ). This is so because Tf^ 3 and T^ 3 contain 
one-photon resonant contributions. Since all the terms 
have different signs around E « 1/2 it is impossible to 
estimate the value of 'jf K3 by just looking at one indi- 
vidual term. 

Fig. 21 shows a three-dimensional plot of r yf K3 as a 
function of E and X with Tio = = O.leV, and -Eio = 
2eV. jf K3 has been plotted for all possible values of E 
and X to show that indeed, jf K3 is zero when X = 
and reaches its maximum value when X = 1. In terms of 
the energy ratios, the function is sharply peaked around 
E w 1/2. 7j'- K ' 3 also gets smaller and approaches zero 
as the values of Tio and increase relative to E\q. In 
order to obtain bigger values of jf K3 , the ratio between 
r 10 (and r 20 ) and E 10 has to be minimized. 

In order to compare r yf K3 and Im(^ res ) we plot the 
two functions together in Fig. [5] along with their ratio. 
r yf K3 and Im("f res ) are both calculated as a function of 
E, with X = 0.5, r 10 = r 20 = O.leV and E 10 = 2.0eV 
(typical values). The typical range of E for many organic 
molecules is also shown. From the plot, it is obvious that 
the approximation that Im(^ res ) dominates jf K3 is not 
a good one for most values of E\ and, in the range of 
typical E, 7/ is smaller than Im(^ res ) by as much as a 
factor of 3. Clearly, there are big differences between the 
two at various values of E. For example, while jf K3 re- 
mains well behaved and goes to zero at E — > 1, /77i(7 res ) 
diverges. Furthermore, while jf K3 = at E = 1/2, 
im(7 res ) w -4xl0~ 43 cm 6 e y- 1 . We conclude that jf K3 
can not be approximated on or near double resonance by 
Im^res). For the typical values of E that are found in 
organic molecules, the two are - on average - different by 
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FIG. 5: A comparison between r yJ K3 and Im(-/ res ) as a 
function of E, with parameters that are typical for organic 
molecules: X = 0.5, E U) = 2eV, V 10 = r 20 = O.leV, and 
N = 1. 

a factor of 3. 

In conclusion, for a given value of E w , in order to max- 
imize jf K3 we must maximize the value of the transi- 
tion moment to the first excited state (X = 1), minimize 
the damping factors (i.e. the inverse radiative lifetime) 
as a fraction of E\q and operate near double resonance, 
where the imaginary part of second hyperpolarizability 
is a maximum (when Eiq sa E20/2). The influence of 
the values of the inverse radiative lifetimes is evaluated 
in Table where we list the maximum values of jf K3 
together with the value of E that maximizes r yf K 3 for dif- 
ferent ratios between the inverse radiative lifetimes and 
Eio- 

III. MAXIMUM VALUE OF 7res FOR 
MOLECULES WITH THREE OR MORE 
ENERGY LEVELS 

As discussed above, it could be argued that the limits 
obtained using a three-level model do not necessarily hold 
for real molecules with more than three energy levels. 
That is, it is in principle possible that when a molecule 
has more levels, the energy distributions allowed by the 
sum-rules result in a higher limit than the values of ^fp ,ax 
listed in table |l| This seems unlikely, since, as we will 
see in the next section, all the measured values of S(lo) 
are much lower than the maximum limit values allowed 
by the sum-rules, indicating that real molecules, with 
more than three available levels, do much more poorly in 
terms of TPA cross-sections than the idealized three-level 
model. 

We argue that the three-level model gives the best 
possible value of <-y™ aa; . I n other words, as the oscilla- 
tor strengths are spread out over more excited levels, the 
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TABLE I: Maximum value of 7 7 as a function of the inverse 
radiative lifetime. E max is the value of E that maximizes 
■yJ K3 . We have used typical values of the parameters for 
organic chromophores of: X = 0.5 and -Bio = 2eV. 
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value of -fY iax decreases. Although this can not be proven 
rigorously for the case of ^J Ko ° using all the SOS con- 
tributing terms, it can be shown that if the values of E\q 
and E20 are known, the explicitly resonant term, 7 res is 
a bounded quantity for any number of levels available 
and that, indeed, the best performance of 7 res is achieved 
when only three states are available. 

To prove our assertion we start with the expression for 
the explicitly resonant term, Eq. which, in the TPA 
regime E20 — 2fauj, can be written as: 



Since all the terms on the left side of the sum are positive, 
it is obvious that fi 2 and /i| will be maximized when: 

^30^30 + £40^40 + ■ • • + E N/ {j, 2 N , Q — 0, (27) 

that is, when only states 0, 1 and 2 are present. 

Now we consider the combination (fc, I) — (1, 1), which 
generates the following sum-rule: 

2 2 2 2 (he) 2 N 

--EW10 + -E31M31 ^ ^-EW'iMjv'i = — ^ — ' 

728) 

which can be rewritten as: 

^21^21+^31^31 + ' ■ ■+EN'lfJ 2 N'i — ( ~ h-EloMlO 

(29) 

All the terms on the left side sum on Eq. [551 are positive 
which implies that will be maximum when: 



^31^31 + ■ • • + E N /ifl Nn — 0, 



(30) 



which again would imply that only the states 0, 1 and 2 
contribute. Therefore, if more states contribute and the 
result is an enhancement of the response, it would nec- 
essarily be due to the contribution of terms exclusive of 
the j res term since we have proven that j res gets smaller 
as the number of levels increases. 

To summarize, we have found that the non-explicitly- 
non-resonant terms can not be neglected in a strict three- 
level system, and that if adding more states results in 
an increase of the TPA response, it will be due to the 
contribution of the non-explicitly-resonant terms in the 
SOS expression. 



IV. APPLICATIONS 
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Since we assume that Eiq and E20 have fixed values, we 
optimize 7 res using /uf and /if 2 as adjustable parameters. 

Now we apply the Thomas-Khun sum-rules to a gen- 
eral system with A^'-levels. The sum-rules can be written 
as: 
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In this section we apply the theoretical results to eval- 
uate the literature values of the measured TPA perfor- 
mance of organic molecules. We do so by comparing 
the measured value of the TPA cross-section of a real 
molecule with the maximum value of the TPA cross- 
section allowed by the sum-rules. The effective number of 
electrons contributing to the nonlinear response is calcu- 
lated by geometrically weighting the number of electrons 
in each conjugated path of the molecule, § 



N, 



eff 



(31) 



where Ski is the Kronecker delta. We generate a sum-rule 
equation by choosing one pair of values of (fc, I). They 
can each be any integer between and N'. 

First we consider the choice of parameters (k, I) — 
(0, 0), which yields, 



E 10 fi 



' E 2 ofJ-lo + Eur 1 ' 



J 30M30 



P „ 2 



2m 



(26) 



where iVj is the number of electrons in the i th conju- 
gated part of the molecule. We note that our method 
for counting electrons is most appropriate for excitations 
that are typical in conjugated molecules. For other sys- 
tems or other types of excitations, an appropriate method 
for counting electrons would need to be used. 

The molecules under analysis are listed in Figures 
{7\ and El together with the effective number of electrons, 



7 



N e ff. Note that the dendrimer structures, shown in Fig- 
ure El from top to bottom shows molecules 20 and 21. 
The structures labelled G-l and G-2 are the i?-groups, 
that when added to molecule 21, makes the dendrimers 
22 and 23. The effective number of electrons are 18.4, 
35.6, 55.8, and 82.7, respectively. 

The TPA cross-section values of molecules 1 to 12 have 
been reported by Rumi and coworkers ;[13j the values for 
molecules 13 to 19 have been reported by Albota and 
coworkers: [l4| and the values for molecules 20 to 23 have 
been reported by Drobizhev and coworkers. |18| The value 
of Eio for each molecule is calculated from the wavelength 
of maximum linear absorption. The values of the inverse 
radiative lifetimes are estimated to be: Tio ~ I^o ~ 
O.leV. 

We note that Frank-Condon considerations suggest 
that a more appropriate estimate for the energy of an ex- 
cited state is the energy at the onset of absorption rather 
than the energy at the absorption peak. For a typical ex- 
cited state energy of 2 eV and width 100 meV, the onset 
is about 50 meV from the peak, leading to a fractional 
error of 2.5 %. For functions that vary as the fifth power 
of the energy, this contributes an error of less than 15 % - 
the low end for the experimental uncertainty of the more 
accurate nonlinear-optical measurements. So, in using 
the peak values to obtain the energies, our calculations 
have an uncertainty of about 15 %; and, the positions of 
the peaks in our dispersion plots could be off by about 
50 meV. 



Influence of the number of electrons 






G-l 



A simple consequence of the Thomas-Khun sum-rules 
is that the two-photon cross-section of a molecule de- 
pends quadratically on the effective number of electrons. 
With the aid of the sum-rules one can rewrite r yf K3 as: 



7 poo 



he 



2 N 2 
E 5 



(32) 



where $ is a dimensionless function that depends on the 
energy ratios, the transition dipole moment ratios and 
the inverse radiative lifetimes but does not depend on 
N e ff. This is a general result that applies to molecules 
with any number of levels. Eq. 1321 together with Eq. Q] 
predict a quadratic dependence on the effective number 
of electrons for the TPA cross-section: 



%) - IN//*- 

^10 



(33) 



Fig. shows a plot of the TPA cross-sections of the 
different molecules as a function of N 2 ^. The data fol- 
low an approximately linear relationship, indicating that 
even though the performance of a molecule is dictated 
by many factors, it correlates with the square of the ef- 
fective electrons. Note that while this scaling appears to 
be approximately universal for the molecules studied, it 
may not hold in other classes of molecules. 




FIG. 6: Molecules measured by Dhrovizhev and coworkers. 
From top to bottom, the molecules are 20 and 21. The 
structures labelled G-l and G-2 are the i?-groups, that when 
added to molecule 21, makes the dendrimers 22 and 23. The 
effective number of electrons are 18.4, 35.6, 55.8, and 82.7, 
respectively [lj| 
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FIG. 7: Molecules measured by Albota and coworkers. [3] FIG. 8: Molecules measured by Rumi and coworkers. [3 
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FIG. 9: S(ui) as a function of iVf// for the collection of mea- 
sured molecules. The data fits to a line with correlation coef- 
ficient of 0.90394. 



B. Predicting the limiting values for real molecules 

Using a three-level model, from the experimental val- 
ues of N e fr, Eio, Tio and T2o, we can find the maximum 
value of 7j K3 allowed by the sum rules, as follows. From 
the theoretical section, the fundamental limit of r yf K3 
occurs when X = 1, so the problem is reduced to find- 
ing the optimal value of E20 that is consistent with the 
sum-rules. We will denote this optimal value as E^ ax . 

The value of E^ ax is substituted into Eq. to obtain 
the maximum value of <5(w), which we will denote S™^ x . 
Since in the resonant TPA regime, -E20 = 27kj, 5™^ x can 
be written as: 



_/ rpmax \ 2 

°3L = r-o-o ill 



hn 2 c 



(34) 



where (-fj 1 ^*) is the orientational average of the fun- 
damental limit of the dressed second hyperpolarizabil- 
ity. In this manner, the maximum values of the TPA 
cross-section are calculated for the different molecules. 
Since most of the molecules are geometrically one- 
dimensional, an isotropic orientational average yields 
(7/) = T//5. The local fields are calculated using the 
Lorentz-Lorenz^^ model with n = 1.4. Fig. 1101 shows 
the two photon performance (the ratio between the ex- 
perimental value of S(u>) and S^ x ). 

All the measured TPA cross-sections fall below the pre- 
dicted limits by an average factor of w 1.5 x 10 -3 . This 
shows that real molecules are far below the maximum 
limit that is allowed by quantum mechanics. Perhaps, 
as in the case of 7/, the low performance ratio is due 
to the existence of many states beyond the dominant 
three-levels, which results in a dilution of the oscillator 
strengths. 

The molecule with highest experimental TPA cross- 
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FIG. 10: TPA performance as a function of E\o for the col- 
lection of molecules studied here. The performance is evalu- 
ated by calculating the ratio between the experimental value 
of 8(10) and S™£ X (maximum value allowed by quantum me- 
chanics in a three-level system) . 



section is 23 with S(uj) — (11000 ± 1100)cm 4 s - a similar 
performance to the rest of the molecules, even though 
its experimental value of S(u>) is 2 orders of magnitude 
higher than any other molecule in the collection. This 
indicates molecule 23 owes its high TPA cross-section 
mostly to the huge number of electrons contributing to 
the nonlinear response, rather than to a better design 
strategy that increases the effectiveness of those elec- 
trons. In fact, in the collection of molecules reported 
by Drobizhev and coworkers, the best performance is 
achieved by molecule 21. Molecule 23 would need to 
have a S(u>) value of approximately 1.4 times bigger than 
its experimentally-determined value in order to be as ef- 
ficient as molecule 21. 

The best performance is achieved by molecule 5. This 
is a somewhat odd result, because molecule 5 belongs to 
the group of homologues 1-5 and the performance of all 
the other homologues is much worse. Rumi and cowork- 
ers report that the experimental uncertainty on molecule 
5 is bigger than for the rest of their reported molecules, 
which could explain the discrepancy. The rest of homo- 
logues 1-4 behave as expected: the values of 5(ui) in- 
crease as the number of electrons is increased, but the 
TPA performance is the same within experimental un- 
certainty and low in comparison to the other molecules. 

This may be better understood if both excited state 
energies are considered. First, the longer the homologue, 
the smaller the value of E10 (red shifting), which in- 
creases the molecular performance due to the l/Ef de- 
pendence. However, for molecules 1 to 4, there is also 
a red shifting of E20 in such a manner that E increases 
with the length of the homologue. The values of E for 
molecules (1, 2, 3, and 4) are respectively (0.802139, 
0.820513, 0.86165, and 0.86165). This increase of E re- 
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suits in a lower performance. However, although E\q is 
red shifted for molecule 5 in comparison with the ho- 
mologues, -E20 stays the same as for molecule 4, result- 
ing in a better-optimized energy ratio (E — 0.812918). 
Therefore, for molecules 1 to 4, the red shifting of E\q 
is accompanied by the red shifting of E20, and the per- 
formance doesn't change appreciably. In contrast, for 
molecule 5, the increase in l/Ef is not counterbalanced 
by an increase in E, so the performance of the molecule 
is increased. 

Molecule 11 which extends the conjugated path of 
molecule 8 shows the same performance factor as its ho- 
mologue and molecule 14 which extends the conjugated 
path of molecule 13 shows a performance factor that is 
the same as its homologue, within experimental uncer- 
tainty. 

Molecules 8 and 13 have the same performance factor 
and also approximately the same value of E\q. Since they 
differ structurally only by the type of radicals attached 
to the central ring, we conclude that the effects of both 
types of radicals are similar in both the TPA performance 
and the wavelength of maximum absorption. Similarly, 
molecules 11 and 14 have the same performance and 
similar values of E\q. Again, they only differ structurally 
by the type of radicals attached to the central ring; so 
we conclude that the effects of both types of radicals are 
similar. 

Molecules 9, 10, 15 and 19 share a similar design strat- 
egy and differ by the type of groups attached to the cen- 
tral rings. This difference has an influence on the TPA 
performance. The best approach appears to be the at- 
tachment of CN radicals, which makes molecule 15 the 
second best molecule in the whole collection. 



C. Comparison with other models 

In the previous section, S™^ x was calculated by find- 
ing the optimal value of E = -E10/-E20 and taking the 
limit X — *■ 1. This requires "tuning" both the values 
of E20 and |/xio| for the given molecule. A more real- 
istic approach is to take the actual experimental values 
of £10 and £20, calculate E, substitute these values into 
S(lu) and calculate the maximum value possible for those 
energy values in the limit X — > 1. We will denote the 
maximum value calculated in this way as S(E) max . 

Thus, for every molecule, we first compute E as the 
ratio between the experimental values of E\q and E20 
and then we calculate the maximum TPA cross-section 
allowed for that energy ratio, S(E) max . This can be com- 
pared with the maximum TPA cross-section possible for 
the same value of E\o (with an optimized value of E), 
&Tl x ■ Figure ITTI plots the ratios of the experimentally- 
determined values normalized to each calculated limiting 
value, S(uj)/5™£ x and S(u) / '5(E) max , as a function of E w 
for each molecule. 

From Fig. El it is clear that the energy ratio, E, is far 
from optimized for all the molecules studied. If it were 



o 

0.01 

CD 



• §(co)/8(E 
" 5(co)/53 L a ' 


^max 




1 

• 


1 

• • 


■ 

• 
• 


1 


- 


■ 


• 


• 






• 


• • 
• 

• 


1 

• 






1 


■ 




■ ■ 


■ 


■ 

li 








■ 


■ 


■ 

■ 


■ 


■ 

■ 


■ 


■ 

■ 

■ 





E io ( eV ) 



FIG. 11: The ratio between the experimental value of 8(u) 
and 8(E) max (maximum TPA cross-section allowed for the 
measured molecular energy ratios) . This is compared with the 
ratio between S(lo) and (maximum TPA cross-section 

allowed for that value of Eio). 



possible to independently "tune" the value of E20 and the 
value of X, the TPA cross-sections could be improved by 
three orders of magnitude. However, if the experimental 
value of E20 stays fixed, the possible improvement for 
typical molecules such as the ones studies here (by tuning 
A) is only an order of magnitude. 

To determine the relevance of the non-explicitly- 
resonant terms, the maximum TPA cross-section allowed 
for a particular value of E for a molecule, 8(E) max , can be 
compared with the maximum value that would be possi- 
ble if the response were determined only by the explicitly 
resonant contribution, j res (Eq. 0). We will denote this 
maximum value S re s- 

We can also calculate the maximum values predicted 
by Kuzyk's approximation technique of extrapolating 
from the off-resonance results. This off-resonance ex- 
trapolation technique (ORET) assumes that jf K3 is 
dominated by j r es, and predicts the following: 



cres 
°ORET 



63.5^ 

vr 



(35) 
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r 2 
1 10 



E20 
T20 



In the off-resonance case, when r„o <C E\$ 
Eq. |23 can be approximated by: 



So RET 




£20/2, 
(36) 
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FIG. 12: The ratio between the experimental value of 8(u) 
and S(E) max , 8 re3 , Soret and S r % ET . 

Fig. E| plots the ratio between the measured TPA 
cross-section, 8(u>), and the various calculated lim- 
its: S(cj)/S res , 8{lo)/5 reTi and 5(u) / 5{E) max , and 
S{to)/S e ^ ETy w hi cn shows that both Soret and S res are 
about an order of magnitude higher than 8(E) max . (Re- 
call that 5{E) max is the fundamental limit of TPA ab- 
sorption using all terms in the three-level model for a 
fixed E that is determined from experiment.) The ORET 
extrapolation using the resonant expression (Eq. I35[l ap- 
pears to be of the same order of magnitude as the ex- 
act calculation while the off-resonant approximation (Eq. 
BrJ)l deviates by about an order of magnitude. Table ITT1 
shows a comparison of the ORET extrapolation and the 
exact results. On average, variations between the ORET 
extrapolation and the exact results are just over 30% - 
the variation being the price of ignoring the non-explicitly 
resonant contributions to "ff K3 - 

Next, we consider extrapolating the off-resonance 
ORET term, given by Eq.lJHl So, while the Soret values 
are consistently larger than S(E) max for each molecule, 
they are approximately related logarithmically by: 

hx(S(E) max ) = XHSoret), (37) 



TABLE II: The values of 8(E) max compared with the results 
of Eq. 1351 for the set of molecules studied here, and the per- 
centage difference between the exact result and the the reso- 
nant ORET formula CEq.lSSI. 
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the limit, which includes the effects of the contributions 
of the non-explicitly-resonant terms to jf K3 - This ap- 
proximation is simpler to use and is more accurate than 
the resonant expression, so provides the researcher with 
a simple expression to evaluate molecules. 



or 

S{E) max = (S OR et)\ (38) 

where A w 0.773. 

To illustrate the validity of Eq. 1381 we list the values 
of 8{E) max and (Soret) ' 773 for the different molecules 
in Table ITTA Within ±20%, Eq. HO is a good approxi- 
mation for most of the molecules. The case of molecule 
19, which has the largest variance, will be discussed later 
in terms of the consistency of the three-model. Thus, in 
our collection of 23 molecules there are only 4 outside the 
±20% range. As such, Eq. [38]is a good approximation of 



D. Consistency of the three- level model 

In the previous section, S(E) max was calculated by 
substituting the experimental values of E w and E 2 q into 
5(ui) and taking the limit X — > 1. Alternatively, with the 
experimental values of E\q 1 E20 and S(u>) we can find the 
value of X that, using the three-level model, would yield 
the same experimental value S(u>). 

So, for a given molecule, we can calculate X 3L , the 
ratio of ImioI/I/^io"! that is derived from the sum-rule- 
restricted sum-over-states three-level model and the ex- 
perimental data. Three different scenarios might result: 
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TABLE III: The values of 5(E) max and (Soret) ' 773 (Eq. 
1381 for the set of molecules studied here, and the percentage 
difference between the exact result and the approximation 
(Eq. EH). 



TABLE IV: The values of E and X' iL for typical TPA 
molecules studied here. E is calculated from the experimen- 
tal values of Eio and E20. X iL is calculated by using it as 
a floating parameter when fitting the TPA cross-section data 
to the three-level model as restricted by the sum-rules. 
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1. The value of X 3L is physically impossible (i.e. an 
imaginary quantity, X 3L > 1, etc). In this case we 
can conclude that the three-level model does not 
describe the molecule. 

2. The value of X 3L is physically allowed but in- 
consistent with the typical values of X for or- 
ganic molecules. For instance, the three-level model 
could predict X 3L = 10 -4 which is possible but 
highly unlikely for real molecules. In this case we 
would also conclude that the three-level model is 
not likely to accurately approximate the molecule's 
TPA cross-section. 

3. The value of X 3L is consistent with physically- 
allowed values of X in a range that is typical for 
organic molecules. In this case, although we can 
not conclude definitively that the three-level model 
accurately describes the molecule, the reasonable- 
ness of the result suggests that it may be suitable 



as a model for real molecules. 



Table IIVI lists the values of X 3L for the collection of 
molecules reviewed here. The predicted values are all 
consistent with what is normally observed for typical 
molecules, except for molecules 16 and 19, which yield 
complex values oi X. Note that these two molecules show 
the largest discrepancies in Table [n] Clearly, for these 
two molecules, the three level model does not do a good 
job of modelling the TPA cross-section. Interestingly, 
these two molecules (excluding 6) have the smallest TPA 
performance ratio, a result that is consistent with our 
assertion that molecules with more contributing excited 
states are less efficient per electron. The rest, however, 
can consistently be described by our simplified three-level 
model. 
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V. CONCLUSIONS 

We have applied the sum rules to the SOS expression 
to calculate the fundamental limits of the dispersion of 
the two-photon absorption cross-section. These results 
apply at all wavelengths, so can be used both on and 
off resonance. Our new rigorous analysis shows that the 
three-level model and truncated sum rules together give 
a reasonable fundamental limit and is a good approxi- 
mation when modelling most molecules. This conclusion 
partly follows from an analysis that applies the sum rules 
to the resonant terms with an infinite number of states, 
and shows that the resonant term is at its maximum in 
the three- level limit. This result is significant because 
it suggests that our ansatz that the three-level model 
accurately describes molecules whose nonlinear suscepti- 
bilities approach the fundamental limit is reasonable. 

We have shown that the non-explicitly- resonant terms, 
which are often ignored, can not be neglected in real 
molecules even near resonance since these terms have a 
profound effect on resonant dispersion. As such, exper- 
iments measuring TPA cross-sections that are analyzed 
near resonance using only the resonant terms may yield 
unreliable conclusions. 

We find that the simple process of extrapolating the 
fundamental limit of the off-resonant two-photon absorp- 
tion cross-section to resonance using only the resonant 
term of the three-level model yields a result that corre- 
lates well with our exactly rigorous results. A logarith- 
mic relationship is proposed that can be used to convert 
the off-resonance extrapolation result to the exact one 
to within an accuracy of about ±20%. This relationship 
is much easier to apply than the full rigorous result and 
should be useful to researchers who are interested in as- 
sessing the TPA performance ratio of a molecule. 

The TPA cross-sections of existing molecules are found 
to be far from the fundamental limit; so, there is room for 



improvement. To reach the fundamental limit would re- 
quire precise control of the energy-level spacing, indepen- 
dently of the transition dipole moments - a task that does 
not appear possible using today's synthetic approaches. 
To make the task within reach, we present an analysis 
of how a molecule can be maximized for a given energy 
ratio, so, researchers only need to focus on adjusting the 
oscillator strength to the first excited state. While this 
approach will not result in molecules with TPA cross- 
sections at the fundamental limit, it can yield a factor 
of ten improvements in the largest TPA cross-sections 
measured. 

Clearly it is best to normalize TPA measurements to 
the fundamental limits when comparing molecules; but, 
we have shown that simply dividing by the square of the 
number of electrons per molecule yields a good metric for 
comparison. Such a simple performance metric is most 
attractive for cases where the transition moments and 
energy levels of a molecule are unknown, which is often 
the case. 

In summary, we report on a mix of theoretical results 
and analysis techniques that will both aide in building 
a deeper understanding of TPA cross-sections and pro- 
vide the experimenter with quick methods for assessing 
measurements. 
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